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Intermolecular forces are known to precipitate adhesion events between sohd bodies. Inspired by 
^\ , a macro-scale experiment showing the hysteretic adhesion of a piece of flexible tape over a plastic 

substrate, we develop here a model of far-field dry adhesion between two flexible sheets interacting 
via a power-law potential. We show that phase transitions from unadhered to adhered states occur 
as dictated by a dimensionless bending parameter representing the ratio of interaction strength to 
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^ ■ bending stiffness. The order of the adhesion transitions, as well as their hysteretic nature, is shown 

to depend on the form of the interaction potential between the flexible sheets. When three or more 
sheets interact, additional geometrical considerations determine the hierarchical or sequential nature 
of the adhesion transitions. 
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INTRODUCTION 



As fabrication technology and nanoscale engineering increase in complexity, it becomes vital to understand small- 
scale interactions between material components. Surface-tension mediated forces play a large role in self-assembly, 
■ not only at the macro-scale [1], but also for micro-electromechanical and nano-electromechanical structures (MEMS 
^ . and NEMS), and as such, a large amount of work has been done in studying the adhesive forces involved [2I, S S]- 
Q \ Carbon nanotubes (CNTs) have attracted significant attention since they were discovered to exhibit novel electrical 
and mechanical properties, and it has been found that CNTs can adhere to each other under the influence of capillary 
forces At these scales, fluid-regulated forces are not the only factors that must be examined. Dispersion (or 

^ \ van der Waals) forces may become more important than at larger scales, and the microscopic intermolecular forces of 
extended media start to have a macroscopic effect on structural stability [sl, I9I]. 

In addition to progress in nanotechnology, many biological systems also display adhesion phenomena whose origins 
can be traced to intermolecular forces. Geckos are known to adhere to smooth surfaces, without any liquid interface. 



The microscopic arrays of hairs, or setae, on the base of the gecko foot are therefore believed to be the source of such 
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■ effective dry adhesion p^, [111, Il2|, [131, [1^ . In cellular biology, cytoskeletal morphogenesis is regulated by complex 
biopolymer networks: Series of long, thin, elastic filaments that form a scaffolding for eukaryotic cells. Mechanical 
properties of macromolecules such as actin filaments or DNA can be measured by force or deflection analysis at small 
j scales, and polymers adsorbed onto a surface or "zipped" to another molecule can be peeled apart by applying optical 



■ tweezers or other external pulling forces 15, 3, 
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Most of the research into interactions between materials at these scales involve close-range, contact, and sometimes 
capillary forces, arid this is the hmit corisidered by mar.y models and experimerits to date Q H However, 
long-range forces due to fixed charge distributions, polar, or even non-polar interactions can lead to adhesion events if 
the right conditions are satisfied. Previous work has characterized the van der Waals attraction between thin flexible 
objects, both theoretically Q, U and experimeritally R. Ir. this paper, we aim to develop ar. understanding 
of the physical mechanisms by which long-range interaction forces compete with elasticity in the adhesion of thin, 
flexible structures. We first introduce and motivate the prototypical system of interest using a macro-scale experiment 
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FIG. 1: (Color online) An example of an adhesion transition between flexible sheets. A piece of adhesive tape is charged 
electrostatically and then moved slowly towards an uncharged surface, with the adhesive side turned away from the surface, 
(a): The charged tape is held far from from the surface and no noticeable bending occurs; (b): Weak bending is exhibited just 
before the critical transition point; (c): At a critical distance the tape moves rapidly towards the surface; (d): Moving the tape 
away from the wall back to its original position shows hysteretic behavior in the shape of the tape. 



showing the hysteretic adhesion of a piece of flexible tape over a plastic substrate. We then develop a model of far- 
field dry adhesion between two elastic, slender sheets interacting via a power-law potential, and study numerically 
their relative adhesion. We uncover that phase transitions from unadhered to adhered states occur as dictated by 
a dimensionless bending parameter representing the ratio of interaction strength to bending stiffness, as well as the 
form of the interaction potential between the flexible sheets. We then generalize our model in order to study the 
interactions between several sheets, and show that additional geometrical considerations determine the hierarchical 
or sequential nature of the adhesion transitions in that case. 



2. MACRO-SCALE EXPERIMENT 



An example of adhesion transition between between elastic bodies due to long-range interactions may be demon- 
strated using everyday materials, namely a piece of adhesive tape and a plastic substrate. As shown in Fig. [H this 
tape can be shown to exhibit complex adhesion properties. The tape is first given a static charge distribution by 
applying it to a piece of plastic and then removing it swiftly. The tape is initially suspended at the distance shown in 
Fig.[T^, sticky side away from an uncharged substrate. As the suspension distance is slightly decreased (Fig.[T]3), the 
tape becomes weakly attracted to the surface. At a critical distance (Fig. [it), the attraction suddenly pulls the tape 
completely to the surface, where it lays flat along the majority of the substrate. As the tape is pulled away from the 
surface, the shape exhibits hysteresis (Fig. [T]i). As the top of the tape is returned to its initial position, the shape 
remains stuck to the surface, even past the distance where it first adhered. 

This simple macro-scale experiment allows us to introduce some qualitative features of the adhesion transition, 
namely a competition between bending and long-range interaction, a sharp transition in shape, and hysteresis. We 
present below a numerical approach to quantify the behavior of similar, but more general, systems. Note that there 
are other characteristics of the macroscopic experiment that we will not attempt to model in this work, namely the 
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FIG. 2: A schematic representation of a two-dimensional cross-section for a system of two flexible sheets (see text for notation). 
Our analysis will be limited to the regime where e = h/ L <^ 1. 

presence of dynamic effects and force due to gravity. 



The system that we study is displayed schematically in Fig. O Two sheets of length L, thickness a <C and 
width d (not shown) are clamped at their left-most edges, separated by a distance and free to interact along their 
lengths. We assume the deformations to be two dimensional, and describe each sheet by the vertical deformation of 
its centerline, denoted y^, with i = 1,2. While DNA and other semiflexible polymers can become kinked, looped, and 
otherwise knotted, this study will be limited to the case where the length ratio e = /i/L <C 1, i.e. the long- wavelength 
limit. 

Under these assumptions the total energy of the system is given by 



where Bi is the bending modulus of the i*^ sheet, and Xi is the horizontal distance. The function V describes 
the interaction potential energy density between the two sheets, as yet unspecified (note that the integration along 
the widths of the sheets has already been performed formally in V). Extremizing this functional yields mechanical 
equilibrium, as shown by the following system of coupled integro-differential equations with boundary terms 



3. 



THEORY 



3.1. Setup 




(1) 




(2) 



(3) 



L L 



y'{'5yi =0 y'{5y[ = 0, 



(4) 



L L 



y'^'Sy2 =0 y'^Sy'^ = 0. 



(5) 



The boundary conditions are set by the physical conditions of the sheets. While there are many possible cases that 
could be examined, we will consider the common physical scenario in which the sheets are fixed and clamped on 
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the left (^i(O) = /i, ^2(0) = 0, ^-(0) = 0, i = 1,2) and the right edge of the sheets are force- and moment-free 
{y'{{L)=y'{'{L)=Q,i = 1,2). 

The potential V can be chosen to describe the physical mechanism responsible for the adhesion between the 
sheets, [sl, [lol, ^^i^ paper, we are considering a general long-range potential of the form V ~ where 

r = yj{xi — X2Y + (^1(2:1) — ^2(^2))^, and n is a positive integer. More specifically, we set 

AnWa- 

where a is the van der Waals-like radius, W is the strength of the interaction, and N is the number of singular modes. 
The sign of determines whether the interaction is attractive or repulsive. We will examine the more specific form 
of this general potential where only two terms remain, an attractive term n = p with Ap = —1 and a repulsive 
term n = q with Aq = +1. This is the familiar Lennard- Jones-like potential that is used to model intermolecular 
interactions [9^. We will also work with the case that a > a, so no "true" contact between the sheets will occur. In 
a related study, Oyharcabal and Frisch [l^ use a van der Waals-like medium range potential with values of p = 3 
and q = 9 to model the attraction between a thin filament and a nonpolar substrate. Other examples include the van 
der Waals interaction between two filaments (p = 6, = 12), polarized attraction between two sheets {p = 2, q > p), 
Coulombic attraction {p = 1^ q > p), and many others (see Ref. \^ for a review). In fact, a surface with an arbitrary 
charge distribution can be represented by a standard multipole expansion, and in a suitable far-field regime a charged 
polymer or conducting elastic sheet can be modeled by this potential as well. Very generally, by specifying the values 
of p and q, any number of potential interactions can be represented, except in the rather exceptional cases in which 
a power-law potential model is insufficient. 

3.2. Dimensionless Equations 

The system described by Eqs. (|2j)-(|5j) is non-dimensionalized by scaling the vertical displacements by /i, and hor- 
izontal distances by L. In what follows variables are understood to be dimensionless. In that case, Eqs. (|2])-(|5]) 
become 

^r(^i) + n^,i/(xi, 2) - n,,i/(xi, ^; q, 2) = 0, (7) 

(^2) - Up^2i{x2.y;p. 1) + n,,2/(^2, i) = o, (8) 

^i(0) = l, y[{0) = y'{{l) = y'{\l) = 0, (9) 
7/2(0) = 0, ^^(0) = ^^'(1) = ^^''(1) = 0, (10) 

where 11^^^ = p(j^L'^~^W/ Bi is a dimensionless quantity, and where we have defined the integral /(x^, y] a, k) as 



^ ^ Jo [(x,-x,)2 + eV]^+' 



with y = yi{xi) - ^2(^2). 



3.3. Asymptotics 

We now take advantage of the long wavelength approximation (e <C 1) to simplify the integrals of the form 
I{xi^ y; a, k). Introducing the substitution xi — X2 = we obtain 

l — xi 

T( OA ^ yi{xi) - y2{xi - ei) 

/ xi,^;a,2) = — - / ^^^-^d(^, 12 
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Expanding to leading order in e, 



r/ OA ^ yi{xi) - y2{xi)di , [ ^ \ nQ\ 
1 Z*"" du / 1 \ 

" e«+i[yi(xi)-y2(a;i)]"y_^(l + u2)f+i +%e"+V ^ ^ 

^ 1 ^ ^^^^ 



"+i[yi(a;i)-y2(a;i)]"r(l + f) ' " Ve"+' 

where we have used u = ^/[yi{xi) — y2{xi)]- The other integrals in Eqs. ©-([HI) are evaluated similarly. Physically, 
Eq. (|15|) expresses the fact that, in the long wavelength limit L ^ h, each sheet see the other one as being locally 
flat, and therefore at leading order the integration along the horizontal direction can be performed first. 



3.4. Identical sheets 

Having derived above the general system of equations for two interacting sheets, we now consider the simplified case 
where the sheets are identical. Defining z{x) as the distance between the sheets, z{x) = yi{x) — y2{x), Eqs. (ff))- (fTUl) 
become 

z{0) = 1, z\0) = 0, (17) 
z"{l)=0, z'"{l) = 0, (18) 

where Q = {2WJpL^aP)/{BhPe), p = {a/hy-PJq/Jp and = r((l + a)/2)/r(l + a/2). Note that the divergent 
behavior 2; = is prohibited thanks to the repulsive part of the potential in Eq. (p!6|) . The dimensionless quantity 
which we refer to as the bending parameter, is a measure of the relative importance of the interaction forces to the 
elastic forces, while is a dimensionless van der Waals-like radius with a numerical prefactor. Hence the symmetric 
system is completely described by the four parameters {Q^ /3^p^q}. 



3.5. Numerics 



The symmetric nonlinear system described by Eqs. (p!6|) - (p!8|) is solved numerically on an adaptive grid to an absolute 
error tolerance 10~^^, along with a standard Newton- Raphson shooting method using MATLAB. A continuation 
scheme in Q and /3 allows efficient computation of nearby systems. This treatment is similar to that followed in 
Ref. Q. 



4. ADHESION TRANSITION 
4.1. Main result 

The main result of this paper is illustrated in Fig. [3l For the particular values p = 3^ q = 9 and /3 = O.O5J9/J3, we 
display the nondimensional bending energy of the sheets Eb — i-e. the sum of the first two terms in Eq. ([T]) — as 
a function of the bending parameter ft. Although here we have chosen values for {p, g', the results are similar for 
other values, with some possible qualitative differences highlighted in the sections below. 

The sudden shape changes, quantified by the bending energy, are reminiscent of the behavior observed experimen- 
tally in ^ For Q <C 10~^, the two sheets are essentially free-standing, as indicated in Fig. [3^. As ft is increased 
the sheets are attracted weakly to one another, resulting in a small shape change (Fig. [8)3). At a critical value of 
the bending energy jumps discontinuously and the sheets abruptly snap together (Fig. [3t). As ^} increases further. 
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n 

FIG. 3: (Color online) Dimensionless bending energy Eb/Eq as a function of the bending parameter, with p = 3, g = 9 and 
P = O.OSJg/Js- Here is a typical bending energy, Eq = Bh^/L^. Representative shapes of the two interacting sheets are 
shown in the different regions. In Fig. [3^, the sheets are essentially straight. Fig. [Sj) shows a slightly bent state due to weak 
attraction between the sheets. Fig.fS]:: Past a critical value of the sheets abruptly adhere to one another. As Q is decreased, 
the sheets retain their adhered character, although the shapes change, as seen in Fig. [Sji (dashed lines indicate shape from 
Fig. [3]:). There is also another sharp transition as Q is decreased even more, and the sheets detach into a bent arc- like shape 
(Fig. [3^). As Q is decreased further still there is a final sharp transition back to the original weakly attracted shapes shown in 
Fig.Et. 



the sheets become more tightly bound, with the undamped portion becoming smaller. If Q is then decreased, the 
system exhibits hysteretic behavior, with the sheets remaining adhered as shown in Fig.[3li (dashed lines indicate the 
previous shape from Fig. [3b). As Q is decreased further still there is another discontinuity in the energy, and the 
sheets once again take on a qualitatively different shape, arc-like, as displayed in Fig. [3^. Remarkably, there exists a 
second (smaller) hysteresis loop on this branch of the energy profile. The jump in the energy at this second hysteresis 
corresponds to a large change in the contact between the sheet end points and the slope of the sheets. Finally, for 
decreasing 1], the sheets return to the positions shown by Fig. [3^ via another sharp transition. 

We also plot in Fig. [H the dimensionless force, ^i'(O), and moment, ^i(O), necessary to apply to the left edge of 
the first sheet to maintain it clamped. For the values of {p^q^/3} considered (6, 12, and O.I5J12/J67 respectively), the 
system models two filaments interacting via van der Waals forces. Much like the shapes themselves, the forces and 
moments undergo sharp transitions and exhibits hysteresis. Note that if the sheets were free to interact they would 
adhere along their entire length, and a force would need to be applied to one end in order to peel them apart. In 
essence the same effect is seen in our system. Hysteresis is known to occur in the strong loading of cantilevers 
and recent experimental investigations into the peeling of CNTs from a substrate have reported results qualitatively 
similar to ours |21[. 
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FIG. 4: (Color online) Dimensionless force, 2/^(0) 5 and moment, (0), on the left side of the first sheet as a function of the 
bending parameter, Q, for values of p = 6, q = 12, and [3 = 0.15Ji2/J6- Note the lack of a second hysteresis loop for these 
values of the parameters. 



4.2. Phase transition and physical quantities 



By viewing the shape-change as a phase transition, where the control parameter is ^ instead of temperature, we 
can borrow several concepts from statistical physics in order to further characterize our model system. The natural 
order parameter to assign is the distance between the sheets, z{x)^ as z{x) = 1 denotes totally unadhered sheets 
and z{x) = (3^1 ^^~p) corresponds to complete adhesion (see Eq. [16]). A natural analogy exists between the energy 
functional given by Eq. ([T]) and a one-dimensional magnetic system with two-component spin subject to an external 
field 0. In our system, there is an energy penalty associated with deforming the sheets (analogously, misaligning 
spins), and there is an interacting field that acts to order the system (analogously, the external magnetic field). It 
is known that even at zero temperature, the magnetic system displays a phase transition at a critical value of the 
ordering field (except in the thermodynamic limit of the sheet length L 00), and as such we could expect such 
behavior from our system as the relative field strength (i.e. Vt) is increased. 

In analogy to the external field of the magnetic system, we define an interaction energy density given by ^{x) = 
Q { — l/z^ + /3/z^)^ and discuss the qualitative changes that govern the phase behavior of the system by studying the 
minima in the free energy (see Ref. [24] for a textbook treatment). When a local minimum appears or disappears 
along the length of the sheets we can expect a change in shape, and whether this change is dramatic or smooth will 
correspond to a first- or second-order phase transition (first-order when dE/dCl is discontinuous, second-order when 
d'^E/dVfi is discontinuous). 

In Fig. [5l we display a representative sampling of the energy densities (bending and interaction energies) with their 
associated shapes, for p = 3, = 9, and j3 = O.O5J9/J3. We see qualitatively different energy densities, confirming 
the transitions between three different phases. In Fig. [5^, the sheets store little elastic energy and are only weakly 
attracted. Past the critical adhesion point, the shape as displayed in Fig.[5)3 now shows a large energetic favorability 
from the interaction force, with large deformation energy penalties on the left edges of the sheets and at the end of the 
un-adhered length. As the bending parameter is increased further, the spatial location of the energy minima shifts, as 
displayed by Fig. [5b. The length of the adhered region (or domain wall), ^, increases with ^7 as (5 ~ 1^"^/^, as expected 
from the boundary layer scaling arising from Eq. ([T6|) . This scaling is confirmed by our numerical simulations (not 
reproduced here). Decreasing ^7, and coming down on the hysteresis loop, at a lower critical value of Vt the adhered 
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FIG. 5: (Color online) Bending energy density, [y" {x)]^ ^ and interaction energy density, ^(x), for various sheet shapes (top) 
along the hysteresis loop for p = 3, g = 9, and jS = O.OSJg / J^'- (a) Weak bending (Q = 0.5), with ends of the sheets approximately 
straight (zero bending energy); (b): Tightly clamped configuration (Q = 0.89) with local minimum in the interaction energy 
density denoting the end of the undamped region. The local minimum in the interaction energy is in the same vicinity, but 
not at the same point; (c): Hysteretic clamped shape {Q. = 0.5), with vertical dashed line indicating the position of the local 
maximum of bending energy density, and the vertical dotted line indicates where this local maximum was in Fig. [Sj). The 
red (outer horizontal dashed) line in the upper inset denotes the weakly bent shape from (a), while the blue (inner horizontal 
dashed) line denotes the shape from (b). (d): Arc-shape (Q = 0.11), with a notable small local maximum in ^(x) near the end 
point indicating localized adhesion. 




sheets become arc-like, and the second local minimum in ^(x) disappears (see Fig. [5]i). 



4.3. General Behavior 

Although our investigation to this point has considered particular values for the parameters g'}, qualitative 

changes in the hysteresis and transition behavior can be obtained for different values of these parameters. Not only 
can the hysteresis region be made to shrink, but it can also disappear entirely. In addition, while the transitions seen 
so far have been first-order, by tuning the model parameters this transition can be made to become second-order. 

These different behaviors are illustrated in Fig. [6l where for convenience we have introduced the parameter h = 
PJp/Jq. As we saw above, there are three characteristic shapes for the sheets that we will denote as weakly bent 
(W), adhered/clamped (C), and arc-shaped (A). The areas that exhibit each of these shapes are depicted in Fig. [H 
In Fig. [6^1 and b, we display the total energy (bending plus attraction) for a system with p = 6, = 12. These values 
model the attraction between two thin non-polar filaments, and the characteristic shapes seen are similar to those in 
Fig.m If h is increased, representing an increase in the minimum adhesion distance between the fibers, the arc-shapes 
and associated phase transition vanish. Similarly for the case where p = 3 and q = ^ pj[, increasing the value of h 
causes the hysteresis region to shrink and the arc-shape vanishes (Figs.[6t and[6]i). Furthermore, the hysteresis region 
can be made to disappear completely for p = 2 and q = A (Fig. [6^ and[6t). 
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FIG. 6: (Color online) A comparison of total energies for the system (bending + attraction energy) for different values of p, 
and (3. For convenience we have introduced the parameter b = f3Jp/Jq. (a): For p = 6 and q = 12, these values model 
the attraction between two thin non-polar filaments. Dotted lines denote sharp phase transitions, while dashed denote smooth 
transitions. This system displays all three characteristic shapes; (b): By changing the value of b so that the sheets do not come 
as close during adhesion, the tightly clamped region is seen to disappear; (c) and (d): For p = 3 and ^ = 9, further qualitative 
changes shrink the hysteresis region and cause the disappearance of characteristic shapes; (e) and (f): In the case where p — 2 
and ^ = 4, the hysteresis can disappear completely even in the presence of a transition. 

The transformation of a first-order transition into a second-order transition indicates that there may be a cusp 
catastrophe in the parameter space we are exploring [25]. For the case p = 3 and g = 9, we plot on Fig. [71 the bending 
energy landscape as both and h = jSJ^,/ J9 are varied. For small values of /3, the first-order nature of the transition 
is apparent. For a fixed value of Vt^ increasing j3 decreases the area of hysteresis, until a critical value is reached where 
the hysteretic behavior vanishes completely. 

4.4. Non- identical sheets 

We have considered so far the case where the two sheets are identical. If we allow instead their bending rigidity to be 
different, we now have two parameters, Vti = {2W JpL^a^) / {BiK'^e) for i = 1,2. In this case, the difference in rigidity 
causes a change in both the equilibrium point of adhesion and the critical values of The different shapes can be 
characterized by a phase diagram that maps the transition points for different values of the two bending parameters, 
as displayed on Fig. [51 Three distinct regions exist: Mult ist ability of both arc-shapes and weak bending (MS), arc 
shapes exclusively (A), and adhered or clamped states (C). As shown in Fig. [8^, for small values of 0.2 and l^i, weak 
attraction occurs, unless the hysteretic regime has been entered, in which case there will be adhesion (not shown). As 
Vti is increased past the fixed value of 0.2^ the sheets adhere in an asymmetric arc shape (Fig. [SJd). This transition is 
first-order, and the sheets snap together. As Vti crosses the second-order transition threshold (dashed line in Fig. [8]), 
the shapes smoothly evolve into a clamped phase (Fig. [8fc). In this final clamped state, the position of the adhered 
portion of the sheets depends on the relative values of the bending parameters (Fig. [8]i). 
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FIG. 7: (Color online) Bending energy landscape for p = 3, = 9 over a range of p and Q. Note that there is a distinct cusp 
in this parameter space, indicative of the "catastrophic" behavior that is associated with first order phase transitions. 



4.5. Adhesion of three sheets 



We now consider the adhesion transition for an array of multiple sheets, and illustrate the complexity and richness of 
the system considered on a few examples. Assuming nearest-sheet interaction for simplicity, with identical potentials, 
we can easily extend the modeling approach offered above to the case of N interacting sheets [26[. The equation of 
shape for the ith sheet, 1 < i < is then given by 



3 



0, 



(19) 



^2 




FIG. 8: (Color online) Phase diagram of sheets with asymmetric bending parameter for p = 3, ^ = 9, and (3 — O.I5J3/J9. There 
is a region of multiple stability (MS) where arc-shapes and weak bending are both possible, arc-shapes only (A) and tightly 
clamped (C) sheets, (a): In this region, the sheets are either bending weakly, or in the hysteretic case they exhibit arc-shapes; 
(b): As Qi is increased, the sheets snap together into an arc-shape (a first-order transition); (c) As Q.i increases further, there 
is a smooth variation from arc-shapes into clamped [i.e. a second-order transition, dashed line); (d) If Qi remains fixed, and 
Q.2 is decreased, the bottom sheet will become relatively more rigid, producing a net shift in the equilibrium adhesion position 
(the dashed lines indicate the adhesion shape in Fig. [5]:). 
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FIG. 9: (Color online) Sequential adhesion transition for three sheets: (1,1,1) (1,2) (3). Two consecutive first-order 
transitions occurring for the total bending energy of the sheets as a function of the (identical) bending parameters in the case 
p = 3, ^ = 9 and P = O.35J9/J3. An inherent asymmetry has been introduced in that the middle sheet is pinned slightly 
closer to the bottom sheet than to the top one (five percent difference in height), (a): Weak attraction; (b) As the bending 
parameters are increased, a first-order adhesion transition takes place where two of the three sheets adhere; (c): A second 
first-order transition occurs when the three sheets adhere. Both transitions display hysteresis. 

where the sum on j runs over nearest neighbors, and fli, p, and P are defined as in the N = 2 case. 

For multiple sheets, any asymmetry in the system now plays a role in determining the order in which sheets adhere 
to one another. For three sheets, unless there is perfect symmetry between the top and bottom sheet, the adhesion 
events always occur in a sequential fashion (see Fig. [9land [TQ|) . Specifically, two sheets first come together, and then 
adhere to the third sheet for a further increase in the relevant bending parameter. By locating the middle beam 
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FIG. 10: (Color online) Same as Fig. (9] but for /3 = OAbJg/Js- In contrast with the case depicted in Fig. [9l the first adhesion 
transition is now of second order. 
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FIG. 11: (Color online) Sequential adhesion transition for four sheets: (1,1,1,1) (1,1,2) (1,3) (4). Total bending 
energy as a function of the (identical) bending parameters for p = 1, ^ = 2, and f3 = O.IJ2/J1. The four sheets form first 
an adhered pair, then a triplet, then all four clamp together. Adhesion in this manner is highly dependent on the asymmetry 
of the array of multiple structures. In this case the top sheet is placed at 2/ = 3/i, the second highest at ^ = 1.8/i, the third 
at 2/ = 0.7h and the lowest beam placed at 2/ = 0. (a): Weakly bent state; (b): Adhesion between the lower two sheets; (c): 
Adhesion between the lower three sheets; (d): Adhesion of all four sheets. Inset: Non-monotonic variation of the bending 
energy in the second-highest sheet. 



slightly closer to one of its neighbors, this sequential can be made to occur preferentially between two previously- 
chosen sheets (by changing the clamping distance, the competition between bending and interaction energy changes, 
and in effect one of the bending parameters gets a boost from the geometric asymmetry). For three identical sheets, 
as the bending parameters are increased, there is a transition from the weakly attracted phase (Fig. [9^1 and [TOk ) to 
a regime where two sheets adhere to each other (Fig. [QJd and [TOb) . As the bending parameters are further increased, 
a subsequent transition occurs where all three sheets come together (Fig. [9b andfTOb). As for the N = 2 case, tuning 
the value of P can change the nature of the first adhesion transition (a^b), from first-order (Fig. [9|) to second-order 
(Fig. [To]). The second adhesion transition, however, remains first-order. 



4.6. Adhesion of four sheets 

When a fourth sheet is added, the adhesion transitions can be made to occur in either a hierarchical or a sequential 
fashion. This is illustrated in Figs. [TTl and[T2l 

We show in Fig. [TT] an example of sequential adhesion for four identical sheets, similar to the one discussed in the 
three-sheet case. The values of p, q and as well as the geometric asymmetry, have been chosen so that there is 
a mix of first- and second-order transitions (see figure captions). The relative distances between the four sheets are 
as follows: The top sheet is pinned at = 3/i, the second-highest at = I.8/1, the third-highest at = 0.7h and 
the final sheet at y = 0. As all four bending parameters are increased at the same rate, the sheets start by a state 
of weak attraction, with more bending exhibited by the sheets that are closer to one another. Past a critical value 
of the bending parameters, there is a second-order transition for these values of the model parameters, and the two 
bottom sheets adhere (Fig. [TTb). As the bending parameters are increased further, another second-order transition 
takes place and the lower three sheets adhere (Fig. [TTb). Finally, a final first-order transition occurs when the four 
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FIG. 12: (Color online) Same as Fig. II II except that the original distances between each sheet is now identical. As a result, the 
adhesion transitions occur in a hierarchical fashion: (1,1,1,1) (2, 2) ^4. 

sheets adhere (Fig. fTTtl). Note the non-monotonic variation of the bending energy in the second- highest sheet (inset 
of Fig.[II]). 

An example of hierarchical adhesion transition is displayed in Fig. [121 where we plot the total bending energy 
profile for a symmetric four-sheet system {i.e. there is no asymmetry in the relative distances between the sheets). As 
the bending parameters are increased, a second-order phase transition leads to adhesion between two pairs of sheets 
(Fig. [T2k). First-order transitions are also possible for other values of the model parameters (not shown here). As the 
bending parameters are increased further, a first-order transition occurs and the four sheets all adhere to one another 
(Fig. [T2b). Remarkably, in this case, the first-order transition does not display any hysteresis. 

5. CONCLUSION 

In this paper we have studied the prototypical dry adhesion problem between flexible sheets or filaments, and focused 
on their morphological transitions. Motivated by a simple macro-scale experiment showing hysteretic adhesion, we 
have introduced a model of dry adhesion between two elastic, slender sheets interacting via a power-law potential, 
and studied numerically the transitions in their conformations. Given a particular form of interaction potential, 
the system is completely described by a single dimensionless parameter quantifying the relative effect of long-range 
attraction and bending rigidity, and governing the nature of the adhesion transitions (first or second-order). We 
have also generalized the model to multiple sheets, showing in particular that additional geometric considerations 
dictate the order in which structure adhere to each other. The physical systems modeled here include the interactions 
between charged sheets, or between nonpolar filaments. Future work will focus on the presence of thermal fluctuations 
allowing the adhered states to "jump" from one state to another. We will also consider the case where the filaments 
are actuated, and will include the effect of hydrodynamic interactions. Finally, using an approach similar to ours, the 
adhesion of three-dimensional structures such as coiled filaments or planar arrays could be investigated. 
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